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In a recent paper ,EP we have reported the ground state 
of a two dimensional system described by a SSH(Su- 
Schrieffer and Heeger)-type Hamiltonlan with a half- 
filled electronic band. The Hamiltoniaro is given by 

H = ~ EK1 - ^(CLaCi+l,^ + h.C.) 
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where the field operators Cij t<T and Cj = a annihilate and 
create an electron with spin a at the site on a 

square lattice subject to the periodic boundary condi- 
tions (PBC) , and the energy is scaled by the transfer in- 
tegral of the regular lattice, the bond variables Xij 's and 
2/ij's in x- and y-directions, respectively, being scaled 
appropriately to involve the electron-lattice coupling co- 
efficient, A representing the dimensionless coupling con- 
stant. In present work we fix the coupling constant at 
A = 0.3, because we are interested only in the properties 
of the ground state which would not depend on the value 
of A. The ground state is accompanied by lattice distor- 
tions which have Fourier components not only with the 
nesting vector Q = (ir/a, itlp>) but also with many other 
wave vectors parallel to Q,t$ 
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(X q , q eW+n a + X 3 * g e- i9 ^>), (4) 
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and a similar expression for y%,j. Here a represents the 
lattice constant. Furthermore, the pattern of these dis- 
tortions is found not unique in the lowest energy state, 
although the structure of the electronic energy spectrum 
and the total energy of the system are the samel? In this 
note we discuss the degrees of degeneracy of the lowest 
energy state for the model described by eq. (1). 

It is clear that there is trivial degeneracies due to the 
symmetry of the system. Three examples of the distor- 
tion patterns in the ground state is shown in terms of the 
wave number dependence of the Fourier components. In 
each graph, the amplitudes of X q ^ q and Y q , q are plotted 
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Fig. 1. The wave number dependence of the Fourier components 
of the lattice distortion. The system size is 32 X 32. For the 
definition of X q>q and Y q , q , see eqs. (2) to (4). 



as functions of q (0 < q < ir /a) on the left-hand side, 
while their arguments are shown on the right-hand side. 
The amplitudes are exactly the same in (a) and (b) of 
Fig. 1 though the arguments look quite different. It is 
not difficult to understand that the difference between 
Figs. 1(a) and 1(b) can be explained by the translation 
symmetry of the system. The behavior of the amplitudes 
in Fig. 1 (c) is completely different from that in Figs. 1(a) 
and 1(b). This fact indicates that there exist non-trivial 
degeneracy in the ground state. From Fig. 1 along with 
many other data which are not shown here, we can sum- 
marize the characteristics of the distortion patterns in 
the ground state as follows; (1) the amplitude of X qA is 
always equal to that of Y q<q , (2) the argument of X qtq 
is equal to that of Y q q or differs from that of Y q , q by tt, 
and (3) the values of and Yk^ are common for all 

the patterns. Thus the variety of the distortion patterns 
in the ground state can be said to be determined by the 
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values of X q> qS with < q < n/a. 

In order to discuss the degrees of degeneracy of the 
ground state, we assume first that the Fourier compo- 
nents of the distortion are restricted to the diagonal ones 
with q_ = q y as numerically confirmed in the previous 
paper,EP and that the displacements y^ a in y-direction 
are automatically determined if the displacements Xjj 
in x-direction are fixed. Then the number of variable 
describing the distortion pattern can be reduced from 
2 x (N - l) 2 to N - 1 for a system with the size N xN. 
It will be allowed to choose x^j (i = 1, 2, • • • , N— 1) with 
a fixed J as a set of N — 1 independent variables. (Note 
the constriction x^j = due to the PBC.) 
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Fig. 2. j for a fixed J is plotted as a function of i. (a), (b) 
and (c) correspond to those of Fig. 1. The value of J is taken 
to be 1 in these cases. Changing J leads to a translation of the 
pattern in the direction of x (or i). 



Figures 2(a) to 2(c) show the distortion x®j for a fixed 
J(= 1) as a function of i, which is constructed from 
X q _qS with < q <.7r /a of the corresponding figures of 
Fig. 1 as in eq. (4) .It It will be clear that Figs. 2(a) and 
2(b) can be completely overlapped by a simple transla- 
tion along i-axis. On the other hand Fig. 2(c) is quite 
different from them. Nevertheless there exists a common 
feature among the three patterns. Namely, (I) xf j van- 
ishes at every second site, and (II) the number of sites at 
which it is positive is the same as the number of sites with 
negative value. Furthermore (III) the absolute values of 
the positive and negative xf/s are equal to each other 
and common to all the three patterns. It will be plausi- 
ble to assume that the absolute value is determined by 
the dimensionless coupling constant A. All the numerical 
data we have obtained show that there is no rule in the 
order of positive and negative xf/s. Thus we can con- 
jecture that the degeneracy of the ground state is given 



by different order of positive and negative xfj's. If N is 
a multiple of four, then xfj will be zero at a half of N 
sites (i = 1, • • • , N) and will be positive at a half of re- 
maining N/2 sites and negative at remaining N/4 sites. 
The degree of non-trivial degeneracy except the trivial 
degeneracy due to arrangement of zero points might be 
said to be N / 2 C N/4 (~ N N / 4 for N > 1). 

In order to check whether the above-mentioned rule, 
(I) to (III), for the distortion pattern is a sufficient con- 
dition for the ground state, we calculate, for the system 
size 32 x 32, the total energy for all the lattice distor- 
tion patterns satisfying the above-mentioned rule, whose 
number is i^Cs = 12, 870. For the amplitude of nonvan- 
ishing xf j, we use the value shown in Fig. 2. In all the 
cases the resulting energy is found equal to the minimum 
energy. Thus the above-mentioned rule can be regarded 
as a sufficient condition. 

It is not easy to show the above rule is a necessary 
condition for the ground state distortion pattern. How- 
ever we can confirm it numerically to some extent. As 
discussed in the previous paper, El we may believe that 
the ground state distortion for the present model has 
only the Fourier components with q x = q y . In this situa- 
tion we can reduce the 2D Schrodinger equation to a ID 
equation in the momentum space. EP In order to perform 
the calculations finding the ground state effectively, we 
adopt this formulation. We start the iterative calcula- 
tion by giving random numbers for the diagonal Fourier 
components of the distortions Xi j 's and yi^ 's and con- 
tinue the iteration up to the 100-th step. The number of 
runs carried out is 100,000. Out of all the runs we could 
obtain 20,445 minimum energy statesjj^ some of them 
are equivalent within the translational shift. We could 
confirm that the above rule is satisfied by all the lowest 
energy states obtained in these calculations. Since we 
find no counter example against the above rule in a num- 
ber of calculations starting from random initial inputs, it 
may be allowed to say that the above-mentioned rule is a 
necessary condition for the ground state distortion pat- 
tern of the present model. We have confirmed the above 
rules for different system sizes, N = 8, 12, 16, 32, 64 and 
128, though the sample number for N = 64 and 128 
are rather restricted, and also for tilted square systems™ 
satisfying PBC's. 

The degeneracy treated in this paper is considered to 
be due to the highly symmetric structure of the model. 
The Peierls distortions in two-dimensional anisotropic 
systems are now under investigation. The results will 
be published elsewhere. 
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